General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



r - 


f 


(NASA-CR- 158236) THE EFFECTS OF CONFINING 
PRESSOBE AMD STBESS DIFFEBEMCE ON STATIC 
FATIGUE GF GRANITE (Laaont-Doher t| 

Geological Observatory) 51 p HC A04/HF AOI 

CSC1 08G G3/46 



M79- 19518 

1 

Onclas 

16412 


The Effects of Confining Pressure 


and Stress Difference 


on Static Fatigue of Granite 


ROBERT L. KRANZ 

Department of Geological Sciences and 
Lamont-Doherty Geological Observatory 
of Columbia University, Palisades, NY 10964 


ABSTRACT 


a 4 n 


Samples of Barre granite were creep tested at room temperature 
at confining pressures up to 2 kilobars. The time to fracture in- 
creased with decreasing stress difference at every pressure, but che 
rate of change of fracture time with respect to the stress difference 
increased with pressure. At 87% of the short-term fracture strength, 
the time to fracture increased from about 4 minutes at atmospheric 
pressure to longer than one day at 2 Kb of pressure. The inelastic 
volumetric strain at the onset of tertiary creep. A, was constant 
within 25% at any particular pressure but increased vith pressure in 
a manner analogous to the increase of strength with pressure. At the 
onset of tertiary creep, the number of cracks and their average 
length increased with pressure. The crack angle and crack length 
spectra were quite similar, however, at each pressure at the onset, of 


tertiary creep. 
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Several theories of static fatigue are all found to adequately 
explain the data. They suggest that the activation enthalpy for the 
stress corrosion process which controls the creep rate increases with 
pressure. The creep strains are best fit by power functions of both 
stress and time. Risking the extrapolation to longer times, the data 
suggest that creep rupture in the upper crust is possible only if de- 
viatoric stresses of several kilobars can be maintained for millions 
of years, if substantial pore pressures can lower the effective 
pressure, or if temperature substantially increases the rate of stress 
corrosion. 
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INTRODUCTION 


The effects of mechanical, chemical and structural variables on 
the fracture strength of rock are well known, having been studied to- 
gether and separately to a great extent. The effects of many of these 
on the time-dependent deformation of rock are less well known. In par- 
ticular, the effect of external or confining pressure on static fatigue 
and creep rupture has not been systematically investigated independent- 
ly of the other variables. 

The problem of the pressure effect on static fatigue is most im- 
portant in the upper crust, becoming subsidiary to temperature effects 
with increasing depth. The transition depth, where creep rupture is no 
longer possible owing to increasing ductility, is unknown. It is sig- 
nificant that dilatancy caused by micro fracturing occurs in crystalline 
rock at pressures up to 8 Kb [Brace et al. , 1966] and temperatures up 
to 400°C [Hadley, 1975a]. Microcracks can coalesce to produce macro- 
fractures in granite at P-T conditions at least up to 5 Kb and 400°C 
[Tullis and Yund, 1377], Unrealistically high differential stresses 
are required at these pressures, however, unless static fatigue, by 
whatever mechanism, has lowered rock strength. 

Early work on this problem by Griggs [1936, 1939], and Robertson 
[1960] on limestone, gypsum and marble concentrated only on transient 
creep where the creep rate persistently decreases. In this stage, 
strain appeared to be proportional to the logarithm of time under load, 
but the proportionality factor was a function of the applied stress 
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difference (maximum stress minus confining pressure). Increasing the 
confining pressure decreased the creep rate while increasing the max- 
imum stress, at any particular pressure. Increased the creep rate. 

The data indicated a change in the mechanism re nonsible for static 
fatigue from dominantly microfracturing below 2 kilobars to some com- 
bination of ductile and brittle deformation at higher pressures. Most 
of these tests were performed with loads that were periodically in- 
cremented so the effect of sustained stress difference was mixed with 
the effects of strain history. 

A large number of uniaxial creep tests (confining pressure equal 
to one atmosphere) have been performed on the silicates, but even now 
few data exist at higher pressures. Most of these are for high tem- 
perature tests where primary interest is centered on steady-state flow; 
for reviews see Carter and Kirby [1978] and Heard [1976]. In this 
regime, the strain rate is nearly constant and is usually described 
by an Arrhenius type equation where the proportionality constant is 
in fact a function of stress. The functional form is usually given as 
the stress difference raised to some power. 

The middle ground, from room temperature and pressure to those 
relevant to the upper crust, is almost barren. Williams and Elizzi 
[1976] have tested gypsum at a constant maximum stress of 400 bars, 
varying confining pressure from 1 to 300 bars. Cogan [1976] has tested 
limestone and shale at low axial loads and pressures up to about 40 
bars. The most complete data come from Wawersik and Brown [1973] and 
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Wawerslk [1972] who tested a granite and sandstone at pressures of 70, 
350 and 690 bars pressure at room temperature. Creep rupture at 
these pressures was observed in only four experiments. Pressure de- 
creased both the primary and secondary strain rates at all stress dif- 
ferences employed. Strain was described as proportional to a power 
function of time under load with the exponent less than unity. 

It is apparent that pressure acts to inhibit creep deformation and 
will increase the time taken to initiate rupture if all other factors 
are constant. Is there a relationship between the pressure effect on 
the short-term strength of rock and its effect on static fatigue? How 
do we account for the pressure effect in existing creep equations and 
theories of static fatigue? With the answers to these questions as its 
goal, an experimental program was initiated several years ago to exam- 
ine both the microscopic and macroscopic effects of pressure and stress 
difference. The current results are given in this paper. 

EXPERIMENTAL PROCEDURE 

By preparing and testing one rock type (in this investigation, 
Barre granite) in a consistent manner, most of the effects of strength 
factors other than pressure and stress difference can be neutralized. 
All samples were cored from two separate blocks in a direction perpen- 
dicular to the rift plane. Block A was the same block used by Kranz 
and Scholz [1977] for their uniaxial experiments; their results will 
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be incorporated within this report. Samples from the two blocks had 
nearly the same pressure dependent fracture strength. Grain size, 
initial crack density and bulk compressibility were similar. Thus 
variability attributable to using separate blocks is discounted. 

Cores were right circular cylinders, 3.46 ± .02 cm in diameter 
ana 8.9 ± .1 cm in length. They were soaked in acetone and dried 
under ambient laboratory conditions. Several days before testing they 
were jacketed with thin copper. Electrical resistance strain gauges 
with effective gauge lengths of both 1.27 cm and 2.54 cm were cemented 
to the jackets without regard to any specific preferred crack orienta- 
tion within the sample. Gauges were aligned to measure strain paral- 
lel (c ) and perpendicular (e ) to the cylinder axis. 
z r 

Volumetric strain (A V/V) was determined in three ways. The 
parallel and two perpendicular gauges were connected directly on the 
jacket; the output was assumed to be proportional to the volumetric 
strain. In addition, the individual outputs from two separate paral- 
lel and perpendicular gauges were externally added electronically to 
measure the volumetric strain. Under stress, rock deforms anisotrop- 
ically, especially while dilating [Hadley, 1975b] and these two sup- 
posedly redundant methods often gave different results. Therefore 
an additional averaging technique was employed. Because the initial- 
ly circular cross section becomes progressively more elliptical, where 
possible the principal strains perpendicular to the cylinder axis were 

determined. Instead of using AV/V - 2s,, + e volumetric strain was 

r z 
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calculated using AV/V ■ e + e + e , where e , and e are the 

r i rz z ri rz 

principal strains of the plane perpendicular to the z axis of the 
cylinder. 

Experimentally, the horizontal 2.54 cm gauge intended to record 
circumferential strain, gave a value somewhere between the two radial 
principal strains, usually near the maximum, so that, volumetric strain 
calculated using it was too high. Volumetric strain from the 1.27 cm 
gauges connected on the sample gave a value close to that calculated 
using the principal strains, probably because the effective gauge 
length of this configuration covers more than half the sample circum- 
ference. Both e and c from the 2.54 cm gauges as well as AV/V from 
r z 

the 1.27 cm gauge configuration are reported. Strains are given with 
respect to the pre-load, hydrostatic state. 

Excepting the uniaxial tests of Kranz and Scholz [1977], all 
samples were deformed in the same hydraulic servo-controlled appara- 
tus. The axial load, applied through a piston to the ends of the sam- 
ple, can be preset to reach and maintain a specified value (within 1 
bar) at a specified loading rate. In addition, a hydrostatic confin- 
ing pressure can be applied and maintained constant independently of 
the axial load. Kerosene was used for the confining medium. The 
loading procedure for each sample, again excepting the uniaxial tests, 
was the same. First, a hydrostatic pressure of 100, 250, 530, 1000 
or 2000 bars was applied to the sample. Next a few additional bars 
were applied to the ends of the sample to seat the piston against it. 
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A load was selected and applied at 100 bars sec -1 . For a fracture 
test, the load was set high enough to achieve fracture. For a creep 
test, the load was set at a value of from 95% to 75% of the stress 
achieved during the fracture test at the corresponding pressure. At 
each pressure, at least one test was made at 87% of the stress differ- 
ence achieved during the fracture test. 

Several tests at 530 bars and 1000 bars of pressure at 87% of the 
corresponding fracture strengths were stopped at, or shortly after, 
the onset of tertiary creep. The sample was unloaded and Internal sec- 
tions of the sample removed for examination with the scanning electron 
microscope (SEM) . These tests were made to augment others by the au- 
thor (Kranz, 1979a] on samples fatigue tested under atmospheric pressure. 
The same sample preparation and crack-counting techniques were used. 

The strains £ , e and A V/V were recorded continuously on strip 
chart recorders. Both e r and e z from the 2.54 cm gauges were recorded 
against the output from an external load cell. Individual values of 
strain are accurate to i 10“ 4 , Load cell output was accurate to 2%. 

No effort was made to increase the load to offset Increases in cross- 
sectional area of the sample which may have reached 2%. Thus actual 
values of stress difference reported are accurate only to within 4%. 
Confining pressure, measured with both a Heise gauge and a BLH pres- 
sure transducer, was maintained constant and known accurately to 
within 1%. 

Times to failure and times to the onset of tertiary creep (where 
strain rate starts increasing) were recorded. In practice, picking 
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the onset of tertiary creep is a natter of laying a straight-edge 
against each of the strain vs. tine records and deciding where the 
slope changes fron concave downward to concave upward. Usually the 
radial strain (e r > curves show this inflection point before either 
the longitudinal strain e or volumetric strain curves. Tertiary creep 
onset tine reported here is the time after loading has ceased to the 
earliest inflection point on any of the strain curves. The accuracy 
of this tine and the failure tine depends on the length of the exper- 
iment because recorder speed was chosen on this basis. Uncertainty 
is estimated at 0.1%. 

OBSERVATIONS AND ANALYSIS 


Table 1 and Table 2 present the data from the fracture and creep 
tests, respectively. Notes at the bottom of each table explain the 
column headings. Fracture strength as a function of confining pres- 
sure is shown in Figure 1. Figure - shows the logarithm of creep rup- 
ture time in seconds as a function of the applied stress difference 
in kilobars at four different pressures. 

Rupture times and theories of static fatigue 

Based upon their experimental work on quartz [Martin, 1972; 
Scholz, 1972] Scholz and Martin proposed that the static fatigue of 
rock could be described by 
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where is the time to failure in seconds, P^q is the partial pres- 
sure of water presumably in the cracks, o is the applied stress in 
kilobars and t Q , a, £ and k are constants. Equation (1) is based on 
data taken at atmospheric pressure and there is no reason to believe 
the constants are independent of pressure or that one can substitute 
stress difference for the a term. The data presented in Table 2 pro- 
vide a valid test of equation (1) , even though the initial P^q was 
not very well controlled; the sample was prepared under ambient 
laboratory conditions. 

Equation (1) can be rewritten in the form 


log (t f ) - A + Bo (2) 

where A ■ In [t P„ _ °] + E/2.303 RT and 

0 »»2v 

B • -k/2.303. Table 3 lists the values of A and B derived by 
fitting equation (2) to the data of Table 2 and also Wawersik f s [1972] 
data for water-saturated Westerly granite (WG) . Uncertainties for the 
slopes, B, are the standard deviations on the fit. Within these 
limits k doesn't change very much. There is a strong pressure depen- 
dence in the A term, however, and this may be due to either a decrease 
in or increase in activation energy. Significantly, the com- 

pletely saturated granite has the lowest A value. 
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Subjected to a hydro* tatlc confining prestur* crack.* close, squeez- 
ing out sons of the air and moisture within the rock. Whet devlatoric 
stresses cause dilatancy, P^ ^ decreases within the cavity spaces which 
open. As there is a snail Increase in the amount of dilatancy at the 
onset of tertiary creep with increasing pressure (A 2 In Table 2), it is 
quite likely that pert of the pressure dependence in the A term of equa- 
tion 2 is a result of lower P^^ values. Equally possible is an in- 
crease in the activation energy term in (1) with pressure. Scholz and 
Martin neglected to include a pressure term In their formulation because 
their tests were all performed at the same pressure. In (1), E should 
be replaced by E* + PV* where E* is the activation energy for the rate 
governing process, P is pressure and V* Is the activation volume. 

Cruden [1974] derived a static fatigue law different in form from 
equation 2. He started by postulating that the strain at the onset of 
tertiary creep is constant. We defer a discussion of this postulate 
until the next section. Further, he assumed that the time to the on- 
set of tertiary creep is proportional to the time to fracture. Anal- 
ysis of the data in Table 2 shows this hypothesis cannot be rejected 
at the 95% confidence level. Thus we assume 

axT - t f (3) 

and using a power law dependence of strain on time and stress. 


e(T,a) - a2 ■ e Q + a30 N T m 
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Combining this with (3) leads to 


t 


f 


*4 


„-*/« 


( 4 ) 


or 


log t. ■ A’ + B' log o 


with A' ■ — log 
m 


r a 1 m (a 2 - e o )~ 

*3 


and B' - — 
m 


(5) 


In his model [Cruden, 1970] N and m can be related to the stress 
corrosion rate and microcrack population. Equation 5 was fit to the 
data in Table 2. The values for A' and B' are given in Table 3 assum- 
ing that a • a i - 03 . 

Does equation 2 or equation 5 fit the data better? Two estimates 
of goodness of fit have been calculated. The first is simply the mean 
square residual (MSR) between the predicted value of log and the 
real value of log t^. The smaller this value is, the "tighter" the 
fit. The second estimator, R 2 , is a measure of how close the model 
comes to fitting the data exactly. If the model were the "correct" 
one, R 2 would be unity. R 2 is the square of the multiple correlation 
coefficient [Draper and S*“‘th, 1966; page 62] and can be used as a 
measure of the adequacy of the model to fit the data. 

Table 3 indicates that bo'.h equation 2 and equation 5 can be made 
to fit the data adequately. vA computer program employing double 
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precision was used.) They are equally "correct" statistically at any 
particular pressure. When all the data (excepting that for WG) are 
used, however, assuming o ■ cj - 03, both equations 2 and 5 are very 
poor because of the additional pressure dependence in the coefficients. 

The similarity of the forms and the fact that stress corrosion 
of cracks is the foundation of both theories, leads or*, to inquire if 
another, more general, model could be as adequate for predictive pur- 
poses. Several models were chosen. The best one had the form 

log t f - a" + B"oi P + C"o 3 S (6) 

This model reduces to equation 2 if 03 < <oj and P is near unity. 

Table 3 indicates that equation 6 is much better than equations 2 or 
5 at fitting all the data but it still has room for Improvement 
(R 2 * .81). Note that P has a value close to one so it is easy to 
see a direct correspondence between the functional forms of the co- 
efficients in equations 6 and 2 . 

Inelastic strain 


It is now well established that creep in brittle rocks is a con- 
sequence of microcracking [Scholz, 1968; Cnden, 1970; Wu and Thomsen, 
1975], even under pressure [Lockner and Byerlee, 1977]. Since crack- 
ing is an inelastic process, attention should be focused on the in- 
elastic part of the recorded strain. Kranz and Scholz [1977] sug- 
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gested that the Inelastic volumetric strain Is a measure of the 
crack density within the rock. 

Tables 1 and 2 list the Inelastic strain components e* and e* 

z r 

as well as the inelastic volumetric strains 4. Figure 3 shows the 
inelastic volumetric strains at each pressure and stress difference. 

For the fracture tests, 4 was recorded near the failure stress, just 
prior to the point where the fidelity of the strain gauge readings 
became questionable; that is, near the point in stress-strain space 
where the instability leading to failure occurred. These points are 
plotted as open symbols in Figure 3 at the a/o c value where they were 
measured. The solid symbols are the 42 values calculated or measured 
at the onset of tertiary creep in each creep experiment, 

Kranz and Scholz [1977] claim that if a critical crack density is 
to be associated with an unstable condition within the rock, the in- 
elastic volumetric strain at the onset of the instability should be 
roughly constant, independent of the applied stress. In Table 4 this 
hypothesis is evaluated at each pressure. The standard deviation of 4 
at each pressure is less than 10% of the average except for the uniax- 
ial data. The correlation coefficient between o/a c and A is small, and 
a T test shows that the hypothesis of a linear correlation between a/o c 
and 4 can be rejected at each pressure. Thus it seems that the insta- 
bility can be marked by a critical dilatant volume, but the magnitude 
of this inelastic volume is a function of pressure, at least up to 2 Kb. 

A comparison of the inelastic volumetric strain near the instab- 
ility in the fracture tests with that at the onset of tertiary creep 
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Is made In Figure 4 as a function of pressure. The value of A In the 
constant-rate fracture tests depends, of course, on the stress level 
chosen to mark the Instability. Even taking this Into consideration 
there Is a rise with pressure In the dilatant volume which the rock 
can stably sustain. Brace [1978, Table 2] and Brace e£ al. [1966, 

Table 3] show a similar result for other rocks, and Wawersik and Brown 
[1973] also found the maximum al 1 owable strain increases as confining 
pressure is raised. Analogously, there is an increase in A with pres- 
sure at the onset of tertiary creep. The dilatant volume at the insta- 
bility is the same within a factor of two in the two different tests. 
Even considering the uncertainty in determining this instability with 
strain gauges, it suggests that a similar critical crack density, or 
some other crack ensemble parameter, must be achieved in both test 
types before the crack ensemble becomes unstable. Whatever happens to 
the crack ensemble in the last 2-5% of a constant-rate test probably 
also occurs in the tertiary phase of static fatigue. 

The process of accumulating crack volume is initially slower in a 
creep test, and apparently allows more strain to build up than in the 
fracture test, Wawersik and Brown [1973] and Scholz and Koczynski 
[1979] have also reported larger accumulated strains in creep tests than 
in faster, constant loading rate tests. In the next section we give ob- 
servational evidence that the increase in dilatant volume sustained with 
pressure is related to size and number of cracks. 

SEM Observations of Cracks 

Ey stepping uniaxial creep experiments at various times and then 
sectioning the sample, Kranz [1979a] was able to study the trend of 
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crack growth as a function of time using the scanning electron micro- 
scope (SEM). Average stress-induced crack lengths increased with time 
as much or more than they did upon loading. In addition, crack Inter- 
action with other cavities seemed to increase with time. It was sug- 
gested that near the onset of tertiary creep crack coalescence was more 
important than individual crack growth in determining rock stability. 
Almost all stress-induced cracks were extensional. 

In this study, the stress-induced crack length and angle spectra 

were compiled at the onset of tertiary creep for samples subjected to 

the same a/cr value (87%) as the uniaxial samples, but at confining 
c 

pressures of 530 and 1000 bars. Analysis was done on block A samples 
only. The data are compared in Figure 5 and tabulated in Table 5. 

Note, in Table 5, that though only half as many traverses across the 
sample area were made on the two samples subjected to pressure, at 
least 30% more cracks were counted. In addition, the average crack 
length in the latter was greater, though crack length was extremely 
variable. The number of cracks with lengths greater than 500 pm was 
greater in the samples subjected to pressure. Average width was only 
slightly greater and also quite variable. The average angle to the max- 
imum stress direction and the spread of crack angles was the same. 

Figure 5 makes these points more clearly. One hundred represen- 
tative cracks are presented. Each unit is one crack counted within 
the appropriate 50 pm * 10° slot. The histograms are quite sim- 
ilar. Most short cracks are within 30° to the maximum stress axis 
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and the longest cracks are within 10°. Those cracks at angles greater 
than 30° are invariably short. At the higher pressures, there are 
more of these. This roughly corresponds to the slight increase in the 
ultimate fracture angle as pressure increases (Table 2). Many of 
these short, off-angle cracks were linking cracks, joining longer ones. 
The histrograms show that at the onset of tertiary creep, the crack 
ensemble has a definitive character and, with the exceptions noted, 
it is the same near the onset of the instability independent of pres- 
sure. 

Creep equations 

Time-dependent deformation, whether in metals, ceramics or rocks, 
has been described using one of three different equation forms: loga- 

rithmic, inverse exponential or a power function with exponent less than 
unity. Often a steady-state factor, proportional to time under load, is 
added on. At least within the rock mechanics literature, there is no 
consensus as to which is the best for predictive purposes. Each has 
its own phenomenological derivation and historical precedent [see re- 
views by Misra and Murrell, 1965 and Cruden, 19711. It therefore 
seemed prudent to try all of them. This is essentially the same ap- 
proach taken by Arfrouz and Harvey [1974]. They concluded that equa- 
tions with time appearing to a power gave the best fit to the data 
for rocks such as coal, limestone and sandstone. Cruden [1971] reached 
the same conclusion. He also concluded that the steady-state term was 


insignificant. 
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All of the above curve-fitting exercises were conducted on uni- 
axial longitudinal strain data. As it seemed possible that strain 
perpendicular to the applied load might not follow the same equation, 
or that the equations might change with pressure, all the strain data, 
at all pressures were fit using least-squares regression techniques. 
Strain was fit as a function of time starting one second after the 
loading and up to the onset of tertiary creep. The results may be 
stated as follows [Kranz, 1979b]: 

1. At all pressures, either 

e(t) - A + Blog t + Ct (7) 

or 

e(t) » A + Bt m + Ct (8) 

could be made to fit the data to a degree much better than 
the accuracy of the data. 

2. Using several criteria for goodness of fit, equation 7 
fit "best" about 50% of the time for the e z component but 
only 25% of the time for either or AV/V for the uni- 
axial experiments. The remainder of the data was best 
fit by equation 8 with or without a Ct term. 

3. Equation 8 without the Ct term almost invariably gave the 
"best fit" at higher pressures for all strains, and when it 
didn't, including the Ct term made it the best. 
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4. In all cases where a term proportional to time was necessary, 
the proportionality constant was of the order of 10" 6 or 
less when time is given in seconds. 

5. The coefficients A, B and C were strongly stress-dependent. 
The stress or pressure dependence of m in (8) remains un- 
clear. 

Since equation 8 was most successful an attempt was made to find 
the and 03 dependence of the coefficient B. Transient creep at 
higher temperatures [Carter and Kirby, 1978] has most often been 
written in the form 


e ■ 8 o o^t m exp (-E/RT) (9) 

Therefore, all the data were refit, using an iterative technique 
[Kranz, 1979b], to the form 

e(t, ai, 03) * a + 8(01 - a3) N t m (10) 

The results for equation 8 are given in Table 6. Equation 10, given in 
Table 7, is only a fair model and there still may be a pressure depen- 
dence within the coefficients. It suggests that 

e(t, 01, 03) * a + 8(ai Nl - a 3 N3 )t m 


(ID 
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might be better. Equation 11 has the additional virtue of providing 
a quick estimate of the relative effects of aj and 93 on the creep 
rate. To determine the best values of a, 3, Ni, N3 and m simultane- 
ously requires 5 pieces of information, or 4 at every time t. As 
0\ and a 3 were not changed during an experiment, equation 11 cannot 
be uniquely defined from this data set. An experiment where 03 is 
incrementally changed, for example, is needed to find the relative 
values of and N 3 . If we accept the m value for each experiment 
as defined by equation 8, however, then it is possible to fit B in 
equation 8 to the form 


B ■ 6(ci Nl - o 3 N3 ) (12) 

using all the data with m as an additional Input. We find Nj * 4.2, 
N 3 * 5 with 8 « -1.1 x 1CT 6 for e.,. and N 3 « .5, N 3 * 1.5 with 
3 ■ 4.15 x 10 ” 3 for e r . The fits are poor (R 2 =.4), however, probably 
because of the large covariance of Nj, N 3 and m. As a result, little 
meaning can be attached to these values, other than to say that the 
strain rate, being proportional to B, is affected more by 03 than by 
ai. 
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Th e effect of pressure 

As this and many other studies have shown, confining pressure 
Increases the differential stress which a rock can sustain before 
fracturing. If a stress difference less than the fracture stress is 
applied, confining pressure increases the fatigue time. If we take 
into account the increase in strength with pressure, is this suffi- 
cient to allow us to predict the static fatigue time? That is, if 
we normalize the applied load to the short-term breaking strength at 
the corresponding pressure, will the data lie on a common fatigue 
curve? 

Figure 6 shows the log of the failure time as a function of the 
applied stress difference normalized to the maximum stress difference 
sustained in the fracture tests at the corresponding pressure. The 
best-fit lines from Figure 2 have been redrawn on this new stress axis. 
Some of the lines intersect at high stress as a result of data scatter 
and different amounts of time to complete loading. The slopes, how- 
ever, indicate that the pressure effect increases with increasing 
pressure. For example, at an applied load of 87% of the fracture 
stress (symbols), the time to failure increases with pressure. Note 
also that the water-saturated granite (WG) is more susceptible to 
stress than the nominally dry granite. 

It might be argued that since microcracking and crack coalescence 
are the basic mechanisms of brittle failure in both the short-term 
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strength teat and tha creep test, then the pressure effect on crack 
growth and coalescence should be the same In both tests. This may 
not be so. Pressure does not affect the nucleatlon of microcracks 
but It does Increase the energy barrier to be overcome for continued 
propagation [Francois and Wllshaw, 1968]. The difference between the 
creep test and rapid fracture test lies In the method of overcoming 
tha energy barrier. In the fracture test, the energy is supplied by 
the continuously augmented deviatoric stress. The cracks are, In 
a sense, driven to greater lengths and interaction. At constant 
stress difference (creep test) cracks propagate into a stable posi- 
tion and stop [Brace and Bombolakis, 1963; Martin, 1972], and furth- 
er growth occurs when the energy barrier is lowered by stress cor- 
rosion reactions at the crack tip [Hillig and Charles, 1965]. Thus 
in the fracture test, crack growth and coalescence are limited pri- 
marily by the stress rate while in the creep test they are limited 
by the corrosion reaction rate or by the rate at which corrosive 
agents can get to the crack tips [Martin, 1972]. 

In order to understand the effect of pressure on static fatigue 
as indicated by Figure 6, microcracking is considered to be a ther- 
mally activated process. Lindholm et al. [1974] show that fracture 
is controlled by such a process. As such, it occurs at a rate pro- 
portional to exp (-AH/RT)f (o) where AH is the activation enthalpy 
associated with the energy barrier at the crack tip, and f(c), a 

N 

function of the applied stress, is usually given as exp(Ka) or a . 
AH - E* + PV* where E* is the activation energy, V* is the activa- 
tion volume and P is pressure. Thus pressure acts to decrease the 
rate at which work is done at the crack tip. 
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One way pressure could affect static fatigue Is through Its effect 
on the stress corrosion process. Stress corrosion Is a thermally act- 
ivated process [Hilllg and Charles, 1965] which occurs at a rate pro- 
portional to 

r /AH' - AH N “[ 

L~ ( ' R T " > ] 

where, to first order, AH’ ■ E* - (a x - a 3 )V* and AH « E q - yV^ 
where y is the crack surface free energy, is the molar volume 
and p is the crack tip radius. Thus devlatoric stress increases 
the corrosion rate while any factor which decreases p or in- 
creases y will decelerate the corrosion process. The effect of pres- 
sure on the corrosion rate is not specifically accounted for in this 
formulation and must come through its effect on AH (by changing p for 
example) or by some additional PV* term in AH’ . Dehart and Liebowitz 
[1968] have shown, for example, that stress corrosion is retarded by 
pressure in some metals, but little other evidence is available. 

Besides increasing the crack propagation energy required and the 
fracture toughness [Schmidt and Huddle, 1977], pressure increases the 
mean normal stress on any plane within the rock. This would be impor- 
tant if friction were a factor in crack growth. The vast majority of 
stress-induced cracks in both constant rate and creep tests are exten- 
sional, not shear cracks, however [Tapponnier and Brace, 1976; Kranz, 
1979a], so the increased normal stress is probably more of a factor in 
closing crack walls. If a crack closes under pressure, the transport 
of fluids and gases through the rock is reduced and this will have an 
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Incidental effect on the corrosion process during creep If new crack 
surfaces cannot be reached by corrosive agents. 

In a creep test, the rate of individual crack growth which is 
controlled by the rate of stress-corrosion, will determine to a large 
extent, the fatigue time. The rate of crack-linking will also affect 
the fatigue time and it is not clear that pressure will have the same 
effect on the crack-linking process as on the individual crack growth 
process. Figure 4 shows that the higher the pressure, the larger the 
crack volume prior to the failure instability. This seems to indicate 
that pressure severely inhibits the crack-linking process. As this 
occurs in both rapid constant-rate tests and creep tests, the rate 
limiting factor for crack-linking cannot be the same as for individual 
crack growth. Apparently, stress-corrosion cracking in the creep 
test or stress-induced cracking in the constant-rate test are initial- 
ly dominant and the crack coalescence rate is only significant near 
the end of an experiment when the critical crack density is approached. 
The mechanism by which pressure inhibits crack linking is unknown but 
may involve a diminishing of the stress concentration around each 
crack as a result of the superposition of the least principal stress, 
or a decrease in the shear stress for the same reason. 

Static fatigue in the upper crust 

To relate the results of this study to the production of faults 
or fractures within the Earth, one must be willing to accept a large 
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extrapolation In tine. If wa assume that fatigue rupture can only 
occur at stresses above the dilatancy point (C' in Table 1, where the 
values are biased toward the low side) then, using equation 6 , we can 
estimate the fatigue time at depth. 

Figure 7 shows the amount of time that dry granite could sustain 
a particular stress state. Qualitatively, the chemical effects of 
pore water and temperature would shift the curves downward. Even con- 
ceding several orders of magnitude to these effects, it is obvious 
that for creep rupture of whole rock to occur in the upper crust, 
large deviatorlc stresses e~e necessary. For example, assuming com- 
pletely saturated granite so that the <73 depth gradient is 166 bars/ 

km, at 3 km depth a stress difference of 2.5 kllobars could be sus- 

15 

tained for over 31 million years (10 sec). At this depth, the tem- 
perature is less than 100°C but Its effect on the stress corrosion 
process could be considerable. Using an activation energy of 1.4 
kcal/mole for granite (Carter and Kirby, 1978; Table 1] the time to 
fracture at 90°C is about 2 orders of magnitude less than that at 
24*C under otherwise identical conditions. 

CONCLUSIONS 


To summarize, the results of this study indicate that the effect 
of pressure is to decrease creep rates, increase fatigue failure time 
and increase the amount of inelastic deformation rock can sustain be- 
fore becoming unstable. The pressure effect can be traced primarily 
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to an increase in the activation enthalpy required for the stress 
corrosion process governing the creep rate- Ancillary pressure ef- 
fects possibly Include the decrease in as dilatancy grows with 
tine and the decrease in crack interaction. The total pressure ef- 
fect on static fatigue is more severe than would be expected from 
a consideration of the pressure effect on fracture strength alone. 

Experiments with conditions closer to those found in the crust 
are needed. In particular, the individual effects of pore pressure 
and temperature under crustal stress regimes need much further study 
in order to characterize the brittle-ductile transition zone. This 
study has shown that large differential stresses can be sustained in 
whole rock at shallow depths for long periods of time but this may 
not be the case at greater depths. 
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TABLE 3. Parameters for Static Fatigue Equations 
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TABLE 4. Evaluation of Inelastic Volumetric Strain Data 
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TABLE 5. SEM Data 
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*20 traverses across the sample were made for the 1 bar samples whereas only 
10 traverses were made across the same area In the other two samples. 
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Figure 2 
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Figure 1, 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


FIGURE CAPTIONS 

Street difference et failure aa a function of confining 
preeaure for faat loading (100 bare sec*' 1 ) fracture 

tests. 

Logarithm of the failure time in aeconda (y axle) vs. 
a treat difference in kllobara (x exit) at confining 
prea auras of (a) 1 bar, (b) 530 bars, (c) , 1000 bars, 
and (d) 2000 bars. Beat fit, least -squares regression 
line la shown. 

Inelastic volumetric strains at stress differences nor- 
malized to the fracture strength at each pressure. Solid 
symbols are 42 at the onset of tertiary creep (from Table 
2), open symbols are A prior to failure (from Table 1). 
Dilatant volumetric strain at the onset of the Instability 
leading to failure as a function of pressure. Open sym- 
bols are for constant-rate fracture tests, closed symbols 
for creep tests. Error bars are standard deviations. 
Histograms of crack length and angle spectra for samples 
subjected to 1 bar, 530 bars and 1000 bars at 872 of the 
corresponding fracture strength. Each unit represents 
one crack within a 5 urn x 30* slot. One hundred repre- 
sentative cracks are shown in each histogram. 
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Figure 6. Log of Che failure time as a function of Che applied load 
normalized Co che short-term breaking scrength at Che 
corresponding pressure. WG Is for water saturated Wester- 
ly granite [Wawersik, 1972]. Because of loading rate dif- 
ferences, to be strictly comparable, the WG curve should 
be shifted downwards about 2 orders of magnitude in time. 
Symbols are the same as in Figure 3. 

Figure 7. Rupture times expected under various maximum and minimum 
stress states. Isochrons are based on extrapolation of 
equation 6 using the parameters listed in Table 3. 
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